Using Lyman-alpha to detect galaxies that leak Lyman continuum

We propose to infer ionising continuum leaking properties of galaxies by looking at their Lyman-alpha line profiles. We carry out Lyman-alpha radiation transfer calculations in two models of HII regions which are porous to ionising continuum escape: 1) the so-called "density bounded" media, in which massive stars produce enough ionising photons to keep the surrounding interstellar medium transparent to the ionising continuum, i.e almost totally ionised, and 2) "riddled ionisation-bounded" media, surrounded by neutral interstellar medium, but with holes, i.e. with a covering factor lower than unity. The Lyman-alpha spectra emergent from these configurations have distinctive features: 1) a "classical" asymmetric redshifted profile in the first case, but with a small shift of the maximum of the profile compare to the systemic redshift (Vpeak < 150 km/s); 2) a main peak at the systemic redshift in the second case (Vpeak = 0 km/s), with, as a consequence, a non-zero Lyman-alpha flux bluewards the systemic redshift. Assuming that in a galaxy leaking ionising photons, the Lyman-alpha component emerging from the leaking star cluster(s) dominates the total Lyman-alpha spectrum, the Lyman-alpha shape may be used as a pre-selection tool to detect Lyman continuum (LyC) leaking galaxies, in objects with well determined systemic redshift, and high spectral resolution Lyman-alpha spectra (R >= 4000). The examination of a sample of 10 local starbursts with high resolution HST-COS Lyman-alpha spectra and known in the literature as LyC leakers or leaking candidates, corroborates our predictions. Observations of Lyman-alpha profiles at high resolution should show definite signatures revealing the escape of Lyman continuum photons from star-forming galaxies.Comment: A&A in pres

Confidence limits of evolutionary synthesis models III. On time-integrated quantities

Evolutionary synthesis models are a fundamental tool to interpret the properties of observed stellar systems. In order to achieve a meaningful comparison between models and real data, it is necessary to calibrate the models themselves, i.e. to evaluate the dispersion due to the discreteness of star formation as well as the possible model errors. In this paper we show that linear interpolations in the log M - log t_k plane, that are customary in the evaluation of isochrones in evolutionary synthesis codes, produce unphysical results. We also show that some of the methods used in the calculation of time-integrated quantities (kinetic energy, and total ejected masses of different elements) may produce unrealistic results. We propose alternative solutions to solve both problems. Moreover, we have quantified the expected dispersion of these quantities due to stochastic effects in stellar populations. As a particular result, we show that the dispersion in the 14N/12C ratio increases with time.Comment: 11 pages, 8 figures, accepted by A&

Insights on star formation histories and physical properties of 1.2≤z≲41.2 \leq z \lesssim 4 Herschel-detected galaxies

We test the impact of using variable star forming histories (SFHs) and the use of the IR luminosity (LIR) as a constrain on the physical parameters of high redshift dusty star-forming galaxies. We explore in particular the stellar properties of galaxies in relation with their location on the SFR-M* diagram. We perform SED fitting of the UV-NIR and FIR emissions of a large sample of GOODS-Herschel galaxies, for which rich multi-wavelength observations are available. We test different SFHs and imposing energy conservation in the SED fitting process, to face issues like the age-extinction degeneracy and produce SEDs consistent with observations. Our models work well for the majority of the sample, with the notable exception of the high LIR end, for which we have indications that our simple energy conservation approach cannot hold true. We find trends in the SFHs fitting our sources depending on stellar mass M* and z. Trends also emerge in the characteristic timescales of the SED models depending on the location on the SFR-M* diagram. We show that whilst using the same available observational data, we can produce galaxies less star-forming than usually inferred, if we allow declining SFHs, while properly reproducing their observables. These sources can be post-starbursts undergoing quenching, and their SFRs are potentially overestimated if inferred from their LIR. Fitting without the IR constrain leads to a strong preference for declining SFHs, while its inclusion increases the preference of rising SFHs, more so at high z, in tentative agreement with the cosmic star formation history. Keeping in mind that the sample is biased towards high LIR, the evolution shaped by our model appears as both bursty (initially) and steady-lasting (later on). The global SFH of the sample follows the cosmic SFH with a small scatter, and is compatible with the "downsizing" scenario of galaxy evolution.Comment: 28 pages, 26 figures, one appendix, Accepted for publication in Astronomy & Astrophysic

On the effective temperature scale of O stars

We rediscuss the temperature of O dwarfs based on new non-LTE line blanketed atmosphere models including stellar winds computed with the CMFGEN code of Hillier & Miller (1998). Compared to the latest calibration of Vacca et al. (1996), the inclusion of line blanketing leads to lower effective temperatures, typically by 4000 to 1500 K for O3 to O9.5 dwarf stars. The dependence of the Teff-scale on stellar and model parameters - such as mass loss, microturbulence, and metallicity - is explored, and model predictions are compared to optical observations of O stars. Even for an SMC metallicity we find a non-negligible effect of line blanketing on the Teff-scale. The temperature reduction implies downward revisions of luminosities by 0.1 dex and Lyman continuum fluxes Q0 by approximately 40% for dwarfs of a given spectral type.Comment: 6 pages, 4 figures. To be published in A&

The BHK Color Diagram: a New Tool to Study Young Stellar Populations

A new method to derive age differences between the various super star clusters observed in starburst galaxies using the two color diagram (B-H) vs (H-K) is presented. This method offers a quick and easy way to differentiate very young and intermediate age stellar populations even if data on extinction are unavailable. In this case, discrimination of regions younger and older than 4 Myr is feasible. With the availability of data on extinction, the time resolution can be improved significantly. The application of the method to the starbursting system Arp 299 is presented. The validity of the method is confirmed by comparing the equivalent width of the H-alpha line with the chronological map of the northern part of NGC 3690.Comment: 32 pages, 7 figures, 1 table, AJ accepte

The Massive Stellar Content in the Starburst NGC3049: A Test for Hot-Star Mode

We have obtained high-spatial resolution ultraviolet and optical STIS spectroscopy and imaging of the metal-rich nuclear starburst in NGC3049. The stellar continuum and the absorption line spectrum in the ultraviolet are used to constrain the massive stellar population. The strong, blueshifted stellar lines of CIV and SiIV detected in the UV spectra indicate a metal-rich, compact, massive (1E6 Msol) cluster of age 3--4 Myr emitting the UV-optical continuum. We find strong evidence against a depletion of massive stars in this metal-rich cluster. The derived age and the upper mass-limit cut-off of the initial mass function are also consistent with the detection of Wolf-Rayet (WR) features at optical wavelengths. As a second independentconstraint on the massive stellar content, the nebular emission-line spectrum is modeled with photoionization codes using stellar spectra from evolutionary synthesis models. However, the nebular lines are badly reproduced by 3--4 Myr instantaneous bursts, as required by the UV line spectrum, when unblanketed WR and/or Kurucz stellar atmospheres are used. The corresponding number of photons above 24 and 54 eV in the synthetic models is too high in comparison with values suggested by the observed line ratios. Since the ionizing spectrum in this regime is dominated by emission from WR stars, this discrepancy between observations and models is most likely the result of incorrect assumptions about the WR stars. Thus we conclude that the nebular spectrum of high-metallicity starbursts is poorly reproduced by models for WR dominated populations. However, the new model set of Smith et al. (2002) with blanketed WR and O atmospheres and adjusted WR temperatures predicts a softer far-UV radiation field, providing a better match to the data.Comment: To be published in ApJ, Dec. issue 17 figures, 3 in gif forma